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We argue that to explain recent resonant tunneling exper- 
iments on crystals of Mni2 and Fes, particularly in the low-T 
limit, one must invoke dynamic nuclear spin and dipolar inter- 
actions. We show the low-T, short-time relaxation will then 
have a \ ftfr form, where r depends on the nuclear T2, on the 
tunneling matrix element A10 between the two lowest levels, 
and on the initial distribution of internal fields in the sample, 
which depends very strongly on sample shape. The results 
are directly applicable to the Teg system. We also give some 
results for the long-time relaxation. 

1. Recent experiments on magnetic relaxation of 
molecular crystals of Mni2 and Fes have found strong 
evidence for quantum tunneling-mediated relaxation at 
low temperatures 0-0] . This evidence comes from strik- 
ing resonances observed in the relaxation rate at certain 
values of applied external magnetic field H (when the 
energy levels of magnetic states on opposite sides of the 
potential barrier match each other). These resonances 
exist both in the low-T limit (when only the two lowest 
electronic energy levels of each molecule are involved), 
and also in the thermally activated regime, when tunnel- 
ing clearly takes place between higher levels. In the Fe 8 
system, the low-T relaxation rate peak at H = is over 
4 orders of magnitude above the rate at H = 0.1T, off 
resonance! 

In these systems the anisotropic potential acting on 
the molecular " giant spins" (of spin quantum number S) 
is dominated by a strong parabolic easy-axis term, of 
form II = — (D/S)Szi tunneling is caused by weak 
transverse perturbations on this. For Fe§, T> ~ 0.27 K, 
whereas in Mni2, 2? ~ 0.61 K\ and S = 10 for both 
molecules 0-0]. Here we concentrate on the relaxation 
at low T (below T c - 2.2 K and - 0.4 K, for Mu 12 
and Fes respectively ||), near H = 0, so that only the 
two lowest levels 1 10} and | — 10) of ^Ti^ are involved 

- they are coupled by a tunneling matrix element A 10, 
which can in principle be calculated Q, but only if all 
relevant transverse couplings are known Jl0[ . Present 
estimates range over several orders of magnitude, but it 
seems unlikely that A10 exceeds 10 -8 K for either system 

- here we will treat it as an independent parameter. 

The low-T, low-H limit is particularly interesting be- 



cause of the following simple reasoning. Recall first that 
a standard phase space argument |Q tells us that the 
phonon-mediated spin relaxation must go like r _1 ~ 
i z [2N{i) + 1], where £ = \S S - £-s\ is the bias between 
the two lowest levels in an applied field, and iV(£) is the 
Bose function. If we ignore the hyperfine and dipolar 
fields, £ = 2gp,gS\H\, and we get t -1 ~ |-ff| 3 ; however 
it is crucial to understand that even if we include these 
fields, this will not change the prediction of a minimum 
in t -1 around H = 0, unless we include their dynamics. 
This is because the typical bias caused by dipolar fields 
alone is £,Dip ~ 0.1 — 0.5 K; with A10 in the nanokelvin 
range the effect of a distribution of static hyperfine and 
dipolar fields in the sample will be simply to force almost 
all molecules off resonance (resonant tunneling requiring 
that £ < A10, in the low-T limit). The only way then 
for the molecules to relax is then via spin-phonon interac- 
tions, and the effect of the spread in static fields is simply 
to smear the minimum around H — 0. Notice that this 
argument holds even if the static fields have a transverse 
component, unless this is ~ 5T (enough to raise A10 
to the experimental resonance linewidth in energy, ie., 
~ 0.1 — IK for the Mni2 and Fes systems fL2|); this enor- 
mous value is ~ 10 2 times the actual hyperfine/dipolar 
field combination! Thus the low-field, low-T experiments 
force us to consider the dynamics of the hyperfine and 
dipolar fields, which, by varying the bias at each molec- 
ular site in time, is capable of continually bringing more 
molecules to resonance. 

At first glance, in the low-T limit only the dynamic 
nuclear fields (ie., hyperfine and nuclear dipolar fields) 
can then play a role in the relaxation- once kT <C T>, all 
dipolar flip-flop processes are frozen out p3| . However we 
shall see that although one needs the rapidly-fluctuating 
hyperfine field to bring molecules initially into resonance, 
after this the gradual adjustment of the dipolar fields 
across the sample, caused by tunneling relaxation, is suf- 
ficient to bring a steady further supply of molecules into 
resonance, and allow continuous relaxation. This pro- 
cess is particularly important in Fes, where the hyper- 
fine couplings are very small. The fluctuating nuclear 
spin field also makes the tunneling incoherent. One can 
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then write down a classical kinetic equation for the mag- 
netisation, whose solution at short time is found below 
to have a square root behaviour, for almost any sample 
shape (although the characteristic time in the decay de- 
pends strongly on the shape). 

We shall thus find that not only do the low-T, low-H 
experimental results force the inclusion of the dynamics 
of the internal fields into the theory - this dynamics also 
leads to a characteristic (and experimentally testable) 
prediction for the form of the relaxation. 

2. We will treat the problem using a Hamiltonian 
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where the first term describes the dipolar-dipolar inter- 
actions between molecules, the second describes tunnel- 
ing, the third couples magnetic molecules to nuclear spins 
{Ik}, and the last term describes interactions between the 
nuclear spins. This is an effective Hamiltonian operating 
in the subspace of the two lowest levels of each molecule; 
we choose the basis set to be \S Z — ±5); t z and t x are 
Pauli matrices, and {i}, {j} label molecular sites. 

We have ignored the transverse part of the dipolar 
coupling, since it can only renormalise Aio in an un- 
measureable way - all flip-flop transitions to states with 
S z 7^ ±5 are frozen out at low T. Nuclear spin effects 
are in principle more subtle - however, since the dipolar 
fields in ([!]) are diagonal and not dynamic unless some 
of the molecules flip, one has only to understand the dy- 
namics of individual molecules, coupled to the nuclear 
bath, assuming that dipolar fields are frozen. This prob- 
lem was solved in Ref. [0, sec. 4. If the nuclear T\ is 
long (which it will be at these temperatures, since it is 
driven by dipolar flip-flop processes), then the hypcrfinc 
bias field acting on a given molecule rapidly fluctuates 
at a rate T 2 _1 , and over an energy scale which is also 
roughly - T 2 -1 . Typically T 2 _1 - 1(T 7 - 1(T 5 A', so that 
we expect Aio <C . Thus at short times we can write 
the bias £j(t), at molecular site j, as £j(t) = £j + 5£j(t), 
where £j results from the sum of the quasi-static dipolar 
and hyperfine fields, with only a small rapidly fluctuating 
component 6£j(t), which nevertheless sweeps over a bias 
range much larger than Aio. One then finds |Q that a 
molecule in quasi-static bias £ relaxes incoherently at a 
rate 
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The parameter £o depends on the average number A of 
nuclear spins which co-flip with S. If A < 1, then £o ~ 



T 2 ; in the opposite limit £ ~ -Mljd^l- For both Fe 8 
and Mni2, A < 1 is most likely fll5fl ; in any case, £o <C 
Ed, where Ed is the total dipolar coupling from nearest 
neighbour molecules, and the exact value of £o will not 
be too important. 

We now define a normalised 1-molecule distribution 
function P a (£, r; i), with J2 a f d U drP a (£, r;i) = 1. It 
gives the probability of finding a molecule at position 
r, with polarisation a — ±1 (ie., in state \S Z = ±<S*)), 
having a bias energy £, at time t. Molecules having 
bias energy £ undergo transitions between \S Z = S) 
and \S Z = —S) at a rate given by (||). Flipping these 
molecules then changes the dipolar fields acting on the 
whole ensemble, bringing more molecules into near (or 
away from) resonance, and leading to a self-consistent 
evolution of P a (£) in time. The general solution of this 
problem requires a kinetic equation for P a (£,r;t). 

3. To derive a kinetic equation for P a {t;, f; t), we again 
assume that dipolar and hyperfine fields are frozen (apart 
from the nuclear Ti fluctuations just discussed), unless 
a molecule flips. All kinetics then come from these flips, 
along with the resulting adjustment of the dipolar field. 
We may then derive a kinetic equation in the usual way 
fL6| , by considering the change in P a in a time St, caused 
by molecular flips, at the rate t^ 1 (£), around the sample. 
This yields 

Pa & r) = -r N x (0 [P a (£, r) - P_ Q (£, r)] 
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where P aa i(£,C i*) i s the two-molecule distribu- 
tion, giving the normalized joint probability of finding a 
molecule at site r, in state \a) and in a bias £, whilst an- 
other is at f', in state \a'), and in a bias P^ is linked 
to higher multi-molecule distributions by a BBGKY-like 
hierarchy of equations |l6| . The first term on the right- 
hand side of ([|) describes the local tunneling relaxation, 
and the second non-local term (analogous to a collision 
integral) comes from the change in the dipolar field at r, 
caused by a molecular flip at f; the dipolar interaction 
Voir) = Ed[1 — 3 cos 2 8]fl /r 3 , where ^o is the volume 
of the unit molecular cell, and J df ' integrates over the 
sample volume. 

We will assume that at t = the sample is fully polar- 
ized; the initial relaxation can then be treated in a dilute 
solution approximation for the fraction (1 — M)/2 <C 1 of 
flipped molecules (where M = / d£ J (df/flo)[P+(£, f) — 
P-U,r)] = fd£j{dr/£l )M(£,r)). The bimolecular 
distribution function P^ factorizes in this limit, ie., 
P (2) (l, 2) = P(1)P(2). The simplest case one may study 
is that of an ellipsoidal sample, where the macroscopic 
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demagnetization field is uniform. The field distribution 
around randomly placed dipoles is well-known JT7[ to be 
a Lorentzian up to a high-energy cutoff defined by Ed 
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r d (t) = ¥/5 E D {i-M(t)) 

E(t)=cE D (l-M(t)) , 
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where c is a sample shape dependent coefficient |18 , and 
E(t) is the time-dependent internal field. Then (H) gives 



M(t) = -M(t)— / die 
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At very short times t < to£o/£x> this gives a linear relax- 
ation M(t) = 1 — 2£/to, unobservable because £o/Ed <C 
1. For i 3> tq^o/Ed one gets 



M(i) = - 
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Since ([?]) itself is only valid when 1— M(t) <C 1, we simply 
write 



M(t) ^l-y/t] 
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The square-root behaviour will be observable experi- 
mentally over a wide time range, since -Ed/Co ^S> 1. Note 
also that T s hort is sample shape dependent even assum- 
ing a homogeneous demagnetisation field. If the sample is 
not ellipsoidal, than the above analysis is not correct be- 
cause the problem becomes essentially inhomogeneous. 
We then return to the kinetic equation (0), and notice 
that if the demagnetisation varies on a length scale much 
greater than the average distance between flipped spins, 
then (^) is simply modified to 
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where r<j(r, t) and E(r,t) are defined in terms of M(f, t) 
analogously to (|J) and(||); the solution is then identical 
to ([lC]) except that T s h or t is modified to 
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where N(Q) = J drY, a P a (£, = 0,f;t = 0) is the zraftaZ 
"density of states" for the dipolar field distribution, in- 
tegrated over the whole sample, at bias £ = 0; typically 
N(0) ~ 1/Eom, where Er> m is the average demagnetiza- 
tion field. 



In order to verify these results, and to see when the 
square-root behaviour breaks down, we have performed 
Monte Carlo (MC) simulations of the relaxation for var- 
ious sample geometries, by the usual procedure - during 
each time interval St <C tq one flips molecules with proba- 
bility 1 — exp{— <5£/t/v(£)} and then recalculates the dipo- 
lar field distribution, now altered by the flipped molecules 
(cf. Fig.l). The system size we can simulate is not really 
macroscopic , but finite size corrections clearly do not 
change the two main predictions coming from eqtns (^)- 
(|l~2|), viz., (i) universality of the square-root relaxation at 
short times and (ii) the characteristic dependence in ( |l~2] ) 
of the relaxation time on sample geometry. Clearly, the 
fastest relaxation will be observed in nearly-ellipsoidal 
samples. 

It will thus be very interesting to check in the low-T 
limit for this square-root relaxation, using different sam- 
ple shapes. Our calculations are most immediately appli- 
cable to the Fes system [Q , where the field distribution is 
almost entirely due to dipolar spread [ pp[ . Confirmation 
of our predictions would then provide strong evidence 
for the dynamic relaxation mechanism discussed here. 
We emphasize that at higher T we do not expect y/t re- 
laxation, since then dipolar flip-flop processes interfere, 
and T\ becomes short pi) ; moreover, the magnetisation 
reversal also proceeds via higher levels, through mixed 
activation/tunneling processes. Coupling to the phonon 
bath is then crucial, which essentially changes the theory. 

Finally, we consider the relaxation when t 3> T s hort- 
This problem is greatly complicated by the development 
of intermolecular correlations in P^ 2 \ etc., so that 

one may no longer factorize them. However one way of 
avoiding this experimentally would be to let the system 
substantially relax at high T, then cool to low T; one 
would then be in the long-time relaxation regime, but 
with initial condition arranged to give a factorizable P*- 2 ' . 
In this limit another analytic solution for the homoge- 
neous (ie., ellipsoidal) case can be found from (jij), when 

M « 1 and P®,foZW) = P„(0^a'(Oi onc finds 
exponential relaxation, at a rate 
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E m axT [l + Kln(E max /-ir^ )] 



(13) 



where k ~ 1 is a numerical coefficient, and E max is the 
spread in dipolar fields in this nearly depolarized limit. 
Details of the derivation will be given in a longer paper 
©■ 
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Figure 1: (a) Monte Carlo (MC) simulations of the 
relaxation in two samples, each made from a cubic lattice 
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of molecules. In (a) the relaxation is shown as a function 
of yjt/ro for (1) a cubic sample of (50) 3 molecules, and 
(2) for a sphere of diameter 50 lattice spacings ((b) shows 
the same relaxation as a function of t/ro). The dashed 
line shows \fi behaviour; the sample relaxation shows 
multimolcculc correlation effects once M(t) < 0.93. 

Figure 2 The density of states iV(£) for the distribu- 
tion of bias fields, integrated over the sample (cf. text) 
is shown for time t = (where finite size effects smear 
the zero energy delta function), and for t = O.Ito, for 
the spherical sample. Energy (and density of states) use 
units where £o = 1 and Ed = 20. The fraction of states 
in the resonant window of width £ around zero energy, 
at t = 0, was 0.79 (sphere), and 0.037 (cube); the ratio 
y/0. 79/0. 036 w 4.6 corresponds fairly well to the ratio 
~ 4 between the straight-line slopes in (a). 
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